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Graphene is a promising material for ultrafast and broadband photodetection. Ear¬ 
lier studies addressed the general operation of graphene-based photo-thermoelectric 
devices, and the switching speed, which is limited by the charge carrier cooling time, 
on the order of picoseconds. However, the generation of the photovoltage could occur 
at a much faster time scale, as it is associated with the carrier heating time. Here, 
we measure the photovoltage generation time and find it to be faster than 50 fem¬ 
toseconds. As a proof-of-principle application of this ultrafast photodetector, we use 
graphene to directly measure, electrically, the pulse duration of a sub-50 femtosecond 
laser pulse. The observation that carrier heating is ultrafast suggests that energy from 
absorbed photons can be efficiently transferred to carrier heat. To study this, we ex¬ 
amine the spectral response and find a constant spectral responsivity between 500 and 
1500 nm. This is consistent with efficient electron heating. These results are promising 
for ultrafast femtosecond and broadband photodetector applications. 


Photovoltage generation through the photo¬ 
thermoelectric (PTE) effect occurs when light is focused 
at the interface of monolayer and bilayer graphene, 
or at the interface between regions of graphene with 
different Fermi energies Ep HHS]. In such graphene 
PTE devices - which operate over a large spectral range 
[ZllHl that extends even into the far-infrared [3] - local 
heating of electrons by absorbed light, in combination 
with a difference in Seebeck coefficients between the 
two regions, gives rise to a photo-thermoelectric voltage 
kpTE = {S 2 — Si){Tei — To). Here Si and S 2 are the 
Seebeck coefficients of regions 1 and 2, respectively, Tei 
is the hot electron temperature after photoexcitation 
and electron heating, and Tq the temperature of the 
electrode heat sinks. The performance of PTE graphene 
devices is intimately connected to the dynamics of the 
photoexcited electrons and holes, which have mainly 
been studied in graphene samples through ultrafast 
optical pump-probe measurements [Mni- As shown in 
Fig. la, the dynamics start with i) photo-excitation and 
electron-hole pair generation; followed by ii) electron 
heating through carrier-carrier scattering, in competition 
with lattice heating, which both take place on a sub-100 
fs time scale; and finally in) electron cooling by thermal 
equilibration with the lattice, which takes place on a pi¬ 
cosecond time scale. The effect of the picosecond cooling 
step in) on the switching speed of graphene devices has 
been studied using time-resolved photovoltage scanning 
experiments with ~200 fs time resolution [T5H2n] . These 
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studies showed that the picosecond electron cooling 
time limits the intrinsic photo-switching rate of these 
devices to a few hundred GHz, because faster switching 
would reduce the switching contrast, as the system 
does not have time to return to the ground state. 
Indeed, GHz switching speeds have been demonstrated 
in graphene-based devices pTH25] . 


However, the most crucial aspects of the photo¬ 
thermoelectric response are captured by the heating dy¬ 
namics, as electron heating corresponds to photovoltage 
generation. Additionally, these dynamics determine the 
ultimate intrinsic carrier heating efficiency and the re¬ 
sulting spectral response. Here, we measure the pho¬ 
tovoltage generation time with an unprecedented time 
resolution, and assess its effect on the heating efficiency 
through spectral responsivity measurements. In an ideal 
photo-thermoelectric detector, all the absorbed photon 
energy is transferred rapidly to electron heat (before en¬ 
ergy is lost through other channels). In this case, dou¬ 
bling the photon energy would lead to doubling of the 
photovoltage (see Fig. lb), and would result in a flat 
spectral responsivity Rpc = Ipc/Pexc, where Tpc is the 
generated photocurrent and Pexc is the excitation power. 
In strong contrast, the spectral response of conventional 
semiconductor-based detectors is not flat at all, since it 
is determined by the band gap: Photons with an en¬ 
ergy below the band gap are not absorbed, i.e. i?pc = 0, 
and the excess energy above the band gap typically does 
not lead to an additional photoresponse, i.e. a decreasing 
Rpc with photon energy |26j . 
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Photovoltage generation on a femtosecond time scale 

In order to capture the time scale of photovoltage 
generation, we perform time-resolved photovoltage 
measurements with the highest time resolution to date 
of ~30 fs. We achieve this by using a broadband 
TirSapphire 85 MHz oscillator (center frequency 800 
nm, bandwidth >100 nm) that creates <20 fs pulses, 
and a pulse shaper that corrects for any dispersion 
(and thus pulse stretching) that the pulses pick up on 
their way from the laser to the device (see Fig. Ic and 
Supp. Info for details). The concept of the experiment 
is illustrated in Fig. Id. We use a pair of ultrashort 
laser pulses, where the two pulses are separated by a 
controllable delay time t and focused onto a single spot 
on the device. The two pulses contain different spectral 
components to suppress coherent artifacts. We then 
record the photocurrent, averaged over a large number 
of pulse pairs, as a function of delay t. After absorption 
of the first pulse, the electron temperature T^i rises 
by ATei, and then starts cooling down to Tq. If the 
second pulse arrives before Tei reaches Tq (i.e. when the 
delay between the two pulses is shorter than the cooling 
time), the electron temperature rises again but by less 
than ATei. This is because the electronic heat capacity 
Cel of a degenerate electron gas increases with electron 
temperature [57] and ATei = dQ/Cei{Tei), with 

Qq the heat in the system before photoexcitation and 
AQ the absorbed power from a laser pulse. Thus, the 
additional amount of generated photovoltage by the 
second pulse is lower than when the pulses would each 
contribute independently. As a result, the photovoltage 
as a function of delay time directly reflects the time 
dynamics of the electron temperature and therefore also 
the PTE induced photovoltage. 

We apply our femtosecond photovoltage sensing tech¬ 
nique to a graphene pn-junction device, that consists of a 
bottom and top dual-gated graphene flake (“dual-gated 
device”, see Supp. Info). A scanning photocurrent 
image is shown in Fig. 2a, where the top gate and 
bottom gate voltages are such that a pn-junction is 
formed at the edge of the top gate. Figure 2b shows 
the photocurrent with the laser focused at this position, 
as a function of back-gate and top-gate voltages. The 
multiple sign reversals indicate that the photovoltage is 
generated through the PTE effect [311]. We show the 
photocurrent generated in this device as a function of 
t in Fig. 2c, which clearly shows a dip around t = 0. 
In Fig. 2d-e, we show the normalized photocurrent dip 
A/pc for two combinations of gate voltages, both in the 
pn-regime. The decays on both sides of t = 0 reflect 
the cooling dynamics with a picosecond time scale, as 
observed in Refs. [T5H5D] . Around t = 0 we notice 
that the photocurrent dip A/pc is remarkably sharp, 
which is only possible when the time resolution of the 
complete system (i.e. the laser pulses and the graphene 
photoresponse) is sufficiently high. Any decrease of 


time resolution, either because of longer pulses or due 
to a slower generation of the photovoltage in graphene, 
would lead to a broadening of the apex of the inverted 
v-shape. 

We use two approaches to quantitatively determine 
the time scale of photovoltage creation. In the first 
approach, we characterize the sharpness of the photocur¬ 
rent dip around t = 0 by taking the time-derivative. In 
the Supp. Info we show that in the case of exponential 
heating dynamics, taking the derivative directly recovers 
the heat dynamics: dlpc/dt{t) « (1 — 
where T^eat and TcooI are the electron heating and 
electron cooling time scales, respectively. Figure 2f 
shows the experimentally obtained dlpc/dt{t) together 
with the heating dynamics using Theat = 80 fs, obtained 
from a best fit. To show that we can indeed resolve time 
scales shorter than the time scale that was accessible 
in earlier time-resolved studies [181120] . we also plot the 
dynamics with Theat = 250 fs. This slower heating time 
clearly does not fit the experimental data. In the second 
approach we develop a model for the photovoltage as a 
function of t, based on heating dynamics induced by the 
pulse-pair excitation and including nonlinear heating 
(see Methods). We show the data for the first (second) 
gate voltage combination in Fig. 2d(e), together with 
the modeled photocurrent change AIpc{t). We find 
excellent agreement between data and model for Theat = 
50 (80) fs and a cooling time of TcooI = 1-3 (1.5) ps. As 
an illustration we show the model for Theat = 250 fs, 
which is clearly in strong disagreement with the data. 
We thus conclude that for the pn-junction configuration, 
photovoltage generation occurs within 50 fs. 

We now put this capability of ultrafast photovoltage 
generation into the perspective of an application. The 
switching speed of graphene optoelectronic devices us¬ 
ing the direct photoresponse is limited by the picosec¬ 
ond cooling time, as shown by earlier reports [MIU, 
and limited to a few hundred GHz. However, we envi¬ 
sion femtosecond photosensing applications by exploiting 
the nonlinear heating response and combining graphene 
PTE devices with time-differential operation. Here, we 
provide one proof-of-principle demonstration: we use the 
graphene photodetector to measure in a direct electrical 
signal the duration of an ultrashort femtosecond laser 
pulse. We compare the derivative photocurrent signal 
dlpc/dt{t) with the optical cross correlation signal that is 
measured by overlapping the two laser pulses in a second 
harmonic generation crystal and monitoring the second 
harmonic signal at 400 nm as a function of delay time 
(orange area in Fig. 2f). The agreement shows that our 
graphene PTE device is capable of measuring the pulse 
width of the laser down to time scales below 50 fs, in a 
direct electrical signal, and without the use of non-linear 
crystals. We note that this technique will work over a 
much broader spectral range (from the ultraviolet |16j to 
the terahertz 0 ) than techniques based on two-photon 
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absorption in silicon photodiodes or frequency conversion 
in nonlinear crystals, and has a similar sensitivity. 


Efficient photo-induced carrier heating 

Having established that carrier heating and PTE 
photovoltage generation occur on a femtosecond time 
scale, we now address how this step affects the energy 
conversion efficiency of graphene PTE devices. The 
main question is if the carrier heating is fast enough 
to out-compete energy loss processes, such as optical 
phonon emission (step ii) in Fig. la). To this end, we 
study the photoresponse for a wide range of photon 
energies, from 0.8 (1500 nm) to 2.5 eV (500 nm). We use 
a laser source (quasi-CW (continuous wave), since the 
pulse duration of 20 ps is larger than the cooling time 
of ~1 ps) with a controllable wavelength and a constant 
excitation power Pexc- In Fig. 3a, we show the (external) 
responsivity i?pc for the dual-gated pn-junction device. 
This spectral response is dominated by the strongly 
wavelength-dependent absorption spectrum of the de¬ 
vice, which is the result of reflections at the oxide-silicon 
interface and the sub-wavelength oxide thickness (see 
Supp. Info). This effect is very similar to the enhanced 
reflection for certain combinations of wavelength and 
oxide thickness, which enhance the contrast of graphene 
in an optical microscope [55]. Indeed, the calculated 
graphene absorption (using Lumerical FDTD Solutions 
software) agrees well with the wavelength-dependent 
responsivity. 

To avoid the strongly wavelength-dependent absorp¬ 
tion, we use a graphene device that is supported by 
a transparent substrate that consists of 1 mm thick 
quartz (“transparent substrate device”, see Supp. Info). 
The flake contains both single- and bilayer graphene, 
with PTE photovoltage generation at the interface 
[5]. Figure 3b shows the responsivity spectrum at the 
monolayer-bilayer interface, together with the measured 
graphene absorption on a similar device. These data are 
obtained from spatially resolved measurements (Fig. 3c), 
which show that the spatial extent of the photoresponse 
does not change with wavelength. The spectral response 
(for constant excitation power) of the device is strikingly 
constant over this broad range of excitation wavelengths. 
The flat Rpc shows that decreasing the number of 
incident photons (by increasing the photon energy) does 
not lead to a decrease in photovoltage: thus a higher 
photon energy gives a larger photovoltage, as in Fig. 
lb. This is in stark contrast with photovoltaic detectors 
based on semiconductors, where the photovoltage gen¬ 
erally decreases for increasing photon energy, meaning 
that excess energy is lost [55]. The flat i?pc at the 
monolayer-bilayer interface is therefore consistent with 
photo-thermoelectric current generation (which also 
applies to the graphene-metal interface; see Supp. Info). 


To understand the flat, broadband i?pc for constant 
absorbed power, we examine what this result means 
for the electron heating ATei = Tgi — To as a function 
of photon energy. First we note that the photovoltage 
scales linearly with power for all wavelengths (see Fig. 
3d), which means that we are in the “weak heating” 
regime where ATei < Tq and the electronic heat ca¬ 
pacity is constant (in contrast to the “strong heating” 
regime in the ultrafast experiment, where the scaling 
is sublinear). The reason for the small heating is that 
we use more than 10 times lower power and quasi-CW 
laser excitation with ~20 ps laser pulses, which is longer 
than the cooling time of ~1 ps (leading to a peak power 
that is 3 orders of magnitude smaller than for ultrafast 
excitation). In this “weak heating” regime, the cooling 
rate is constant |20] . which means that the conversion 
efficiency is not affected by the life time of the hot 
electrons. Furthermore, the Seebeck factor S 2 — Si does 
not change with excitation wavelength. Therefore, from 
the flat i?pc we conclude that the light-induced increase 
in electron temperature ATei at constant power is the 
same for all photon energies. 

This result enables us to assess the efficiency of the 
electron heating (in the “weak heating” regime). We 
examine two alternative ultrafast energy relaxation 
pathways for photo-excited electrons and holes: carrier- 
carrier scattering and optical phonon emission (see Fig. 
la). Graphene optical phonons have an energy of ~0.2 
eV and therefore photoexcited carriers above 0.2 eV 
can in principle relax by emitting a phonon. The faster 
process of these two competing processes will dominate 
the ultrafast energy relaxation. We have determined 
the time scale of carrier-carrier scattering (through the 
photovoltage generation in the “dual-gated device”) 
to be <50 fs. The time scale of phonon emission is 
typically <150 fs [mill], and therefore this does not 
give a definite answer about which ultrafast relaxation 
process dominates. However, from the measured spectral 
responsivity we can extract the heating efficiency. 

We illustrate this by considering two contrasting cases 
(see Fig. 4a): i) dominant coupling of photo-excited 
electrons to optical phonons, with a small fraction of 
the absorbed photon energy converted into hot electrons 
and ii) dominant carrier-carrier scattering |29[ 130] 
with a large fraction of the energy converted into hot 
electrons. In case i, the energy loss rate dE/dt through 
optical phonon emission increases linearly with initial 
electron energy Ei, since it is governed by a constant 
electron-phonon coupling, and an energy-momentum 
scattering phase space that increases linearly with 
energy [T5] (see Fig. 4a). Thus in case i, a larger Ei 
leads to more energy loss to phonons. On the other 
hand, in case ii, the electron temperature scales linearly 
with Ei, because the energy of the photoexcited electron 
is fully transferred to the electron gas. Thus, the role of 
optical phonon emission can be measured by studying 
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the scaling of ATei with Ei. This relationship, extracted 
from the photovoltage measurements, is shown in Fig. 
4b, where we plot the internal quantum efficiency (IQE 
= ). The IQE represents the generated photovolt¬ 

age, normalized by the number of absorbed photons. 
We find (nearly perfect) linear scaling of the IQE with 
Ei, as the linear fits go nearly trough the origin. This 
means that a higher photon energy corresponds to a 
larger photovoltage and thus to a larger ATei, which is 
consistent with terahertz photoconductivity measure¬ 
ments O [H |30], where a terahertz probe provided a 
measurement of the carrier temperature. We therefore 
conclude that the generated photovoltage comes from 
ultrafast, efficient photon-to-electron-heat conversion 
and an unmeasurably small loss to optical phonons. 

To show that carrier heating is indeed efficient, we cal¬ 
culate the expected temperature rise ATei and use the 
measured photovoltage to determine the Seebeck factor 
S 2 —S 1 , which we then compare to its expected value. De¬ 
tails are given in the Methods section. To calculate ATgi 
we use a simple heat equation and assume fully efficient 
carrier heating to find ~0.17 K (Pgxc = 50 fiW, Tq = 300 
K, Ep = 250 meV, Tgooi = I ps). We then use the mea¬ 
sured photovoltage of Fpte = IpcR ~ 10 where R is 
the device resistance (2 kfl). We conclude that S 2 — S 1 « 
65 /iV/K. This is very close to the maximum expected 
value using a charge puddle width of A = 80 meV, which 
gives 90 fiY /K. Having confirmed that carrier heating is 
efficient, we aim to gain insight into the other factors that 
determine the overall energy conversion efficiency of the 
device. One important factor is the hot carrier lifetime, 
which should be long (low cooling rate) to lead to a larger 
photovoltage. We verify this by changing the ambient 
temperature. First we demonstrate that the efficiency of 
electron heating is independent of lattice temperature, 
as we obtain the same linear scaling through the origin 
for a lattice temperature of 40, 100 and 300 K. We then 
note that the overall photovoltage is larger for lower lat¬ 
tice temperatures. This is caused by the longer cooling 
time Tcooi at low temperatures, due to a lower coupling 
between electrons and acoustic phonons [alEOlEI]. The 
Seebeck coefficient decreases with temperature, meaning 
that the overall generated photovoltage is about a factor 
two larger [20] . 


Conclusion 

The unique femtosecond time resolution and the re¬ 
lated high carrier heating efficiency are very encouraging 
results for bias-free (passive) PTE photodetectors. To 
improve the internal quantum efficiency of ~1% (Fig. 
4b), an interesting approach (see Methods for details) 
would be to use ultraclean, defect-free graphene, such 
as in Ref. [32], which enables detector operation with 
a higher Seebeck coefficient, because the electron-hole 
puddle density is lower (see Methods). The small puddle 


width could increase the Seebeck factor S 2 — to ~ 
300 fj,V/K. Measuring at Ep = 50 meV instead of 250 
meV would furthermore lead to a larger ATei, since 
this scales with 1/Ep (due to the lower heat capacity 
upon approaching the Dirac point). With 50 ^W 
excitation a ATei of almost 1 K is feasible (assuming 
an unmodified carrier cooling rate), giving Vpte ~ 300 
/rV. A resistance of i? = I kD [32] would then give 
an IQE of 100% or more. Our results thus show that 
graphene PTE devices exhibit ultrafast, efficient and 
broadband photodetection. Future photodetection and 
light harvesting devices could exploit these graphene 
properties and combine them with the advantageous 
properties of other 2D materials. Future work could 
investigate the effect of the Fermi energy on the heating 
time, the time-resolved signal generation in devices 
based on the bolometric effect [33l |34] , and the effect of 
the substrate on the heat dynamics. 


Methods 

Simulation of the photocurrent dynamics 

Here we describe how we extract the hot electron 
dynamics from the photocurrent dip A/pc(t). This 
analysis draws on the procedure described in Ref. 
[20] . We start with a laser-induced change in electron 
temperature for single pulse excitation assuming a linear 
response: ATeipp(t') = (1 — Here 

t' is the ’real’ time, Theat the electron heating time and 
Tcooi the electron cooling time. We then introduce a 
nonlinearity (any type of nonlinearity works) of the 
form ATgj = -^1-1- ATgipp/Tg — 1 (with Tq the lattice 
temperature) and calculate the integrated photovoltage 
generated by two pulses well separated in time (adding 
up independently): 2 • Vip = 2 • AT^j ipdt'. We 
now follow the same procedure for two-pulse excitation 
where the pulses arrive with a separation of delay 
time t and obtain ATei^ 2 p (linear response) and AT^j 2 p 
(with nonlinearity). This temperature change AT^j 2 p 
is shown in Fig. Id of the main text in ’real’ time for 
two different delay times t. We use these electron tem¬ 
perature dynamics to calculate the integrated voltage 
V 2 p, and repeat this for a range of delay times t. The 
photocurrent dip as a function of delay time is then 
given by AJpc(t) oc 2 • Vip — V 2 p{t). This dip grows 
upon approaching t = 0, but flattens when t < Theat- 
The reason for this is that when the two pulses arrive 
at exactly the same time {t = 0), there is no heating 
yet due to any of the two pulses and they independently 
start leading to electron heating and thus contributing 
to the photovoltage. As a result, the photocurrent dip 
flattens around t = 0. We note that the dip should 
not disappear, since we measure the time-integrated 
photocurrent and once the electrons start heating up, 
the heating induced by the two pulses is no longer 
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independent. Finally, we point out that all dynamics 
are determined by the local (at the junction) carrier 
heating and PTE voltage creation. All dynamics that 
happen after this (propagation of the potential to the 
contacts, signals traveling through the cables, etc.) do 
no affect the results of our specific experiment. To 
measure transit times, one can use two-pulse excitation 
with two different focus spots [55] . 


Calculation of AT, Vpte and S 2 — Si 

We calculate the theoretical and experimental photo¬ 
voltage VpTE in the case of excitation with Pexc = 50 
/iW of light (of which ^2 % is absorbed, See Fig. 3b), in 
the case of the “transparent substrate device”, where the 
photovoltage is created at the interface of monolayer and 
bilayer graphene, at room temperature. We calculate the 
temperature increase from a simple heat equation (in the 
case of a cooling length that is smaller than the spot size) 
to be |6]: 


AT = Tel - To = 


AQ 

2(Py/4:Tr \n 2 • Fcooi 


0.17K . 


( 1 ) 


Here, the laser focus spot size is d « 1.5 ^m (see Fig. 
3c), and the cooling rate is given by Fcooi = where 
a = 3 ^ 2^2 kg, with h, vp and kp the reduced Planck 
constant, the Fermi velocity and Boltzmann’s constant, 
respectively. We use Ep = 0.25 eV and a cooling time 
of 1 ps, which we obtained from ultrafast photocurrent 
measurements on the same device. This yields a cooling 
rate of 0.4 MW/m^K, close to theoretical estimates that 
produce 0.5-5 MW/m^K |5]. To obtain the photovoltage 
we use jS): 


fpTE = IpcR ~ 10 , (2) 


where we use a photocurrent of 5 nA for Texc = 
50 fiW and a device resistance of i? = 2 kH. Using 
HpTE = ATei(S '2 — A^i), we thus obtain (S '2 — Si) = 
65 /K. The maximum Seebeck factor is given by 

{S 2 — S'i)max = |1|, which gives 90 /rV/K at room 

temperature for a charge puddle width of A = 80 meV. 
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Figures 



FIG. 1: Hot electron dynamics and their experimental extraction, a) Schematic representation of the electron 
dynamics in graphene after photoexcitation with two different photon energies (green and blue arrow). After photoexcitation 
and electron-hole pair generation, electron heating and lattice heating take place. The former occurs through carrier-carrier 
scattering and leads to hot electrons [Hi ll.li 114111611201130] , which drive a photo-thermoelectric current Hi]. The latter occurs 
through phonon emission and leads to the generation of a much smaller photovoltage, since the heat capacity of the phonon 
bath is much larger than that of the electron bath 0- b) The hot electrons lead to a local photovoltage through the photo¬ 
thermoelectric effect, with dynamics governed by the hot electron dynamics. If electron heating is efficient, a higher photon 
energy leads to a larger photovoltage (green line), compared to a lower photon energy (blue line), c) Schematic representation 
of the ultrafast photovoltage setup and the “dual-gated graphene device” (more details in the Supp. Info). The setup produces 
two ultrashort pulses, separated in time by a controllable delay time t. A pulse shaper compensates for pulse stretching in 
the objective. The compressed pulses are focused onto the device, which contains a back gate (BG) and a top gate (TG) for 
creation of a pn-junction at the interface between graphene regions of opposite doping. We read out the photocurrent through 
the source (S) and drain (D) contacts, using a pre-amplifier and a lock-in amplifier, synchronised with the optical chopper, 
d) The development of the electron temperature after excitation with an ultrafast pulse pair with t = 1 ps in the situation 
of independent heating due to the two pulses (black area) and heating where the heating due to the second pulse depends 
on the heating due to the first pulse (purple area). The measured photocurrent /pc is proportional to the purple area; the 
photocurrent dip A/pc is proportional to the difference between the black and the purple area. If t = 5 ps the black and purple 
areas are very similar, i.e. A7pc ~ 0. Thus, the photocurrent dip as a function of t reflects the heating dynamics. 
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FIG. 2: Femtosecond sensing of hot electrons, a) Scanning photocurrent image (green-blue color scale) of the “dual¬ 
gated device”, with indications of the edges of the two contacts and the topgate, and the location of the graphene flake (grey 
rectangle). At the edge of the topgate, there is an interface of two graphene regions whose Fermi energy is separately controlled 
by the voltages on the back gate and the top gate, here tuned to the pn-junction regime, b) The photocurrent as a function 
of gate voltages at the position of the yellow dot in panel a, showing clear multiple sign reversals that are indicative of PTE 
current generation HH]. c) The photovoltage as a function of delay time t shows a dip A/pc when the two pulses overlap and 
recovers with dynamics that correpond to the hot electron dynamics. We show the dynamics of the (normalized) time-resolved 
photovcurrent dip A/pc around f = 0 (blue circles) for gate configurations 1 (d) and 2 (e), both in the pn-junction regime. 
The black solid lines describe the model results (see Methods), using a heating time of 50 (80) fs and a cooling time of 1.3 
(1.5) ps for gate conhguration 1 (2). The red lines show the modelled dip with a slower heating time scale of 250 fs, which is 
incompatible with the data. The insets show the data and model results over a larger time range, f) The time derivative of 
the photocurrent dip (blue circles) in gate configuration 2, together with the derivative of the modelled photocurrent dip with 
heating time scales of 80 (black dashed line) and 250 fs (red dashed line). The orange area represents the optically measured 
cross correlate using a nonlinear crystal. 
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FIG. 3: Spectral response, a) Responsivity as a function of excitation wavelength (blue data points, left axis) and modelled 
absorption spectrum (red line, right axis) for the “dual-gated device” in the pn-configuration, measured with a fixed power. 
The inset shows the device layout, with a 300 nm Si02 substrate on top of Si. Error bars are calculated from independent 
measurement scans and represent the 68% confidence interval, b) Responsivity as a fnnction of excitation wavelength (blue 
data points, left axis) and modelled absorption spectrum (red line, right axis) for the “transparent substrate device” at the 
monolayer-bilayer graphene interface (~50 /rW power). Error bars are calculated from independent measurement scans and 
represent the 68% confidence interval. The left inset shows a scanning photovoltage image with photovoltage generation at 
the monolayer-bilayer (SLG-BLG) interface. The right inset shows the device layout, with a 1 mm thick Si02 substrate. The 
agreement between the responsivity and absorption curve shows that a lower number of absorbed photons (shorter wavelength, 
higher photon energy) does not lead to a lower responsivity, consistent with PTE current generation [3]. c) Photocurrent as a 
function of laser spot position for wavelengths between 500 and 1500 nm (offset for clarity), while scanning the laser through 
the interface between monolayer and bilayer graphene, showing constant photovoltage amplitudes and spatial extent. The laser 
focus spot size is 1.5 ± 0.15 fim. d) Power dependence of the photovoltage with the laser focused at the SLG-BLG interface 
(see inset in panel b) for a range of wavelengths (see inset for the wavelength corresponding to each color), showing linear 
scaling, which corresponds to the “weak heating” regime, where ATei < To- The inset shows the fitted power exponent, which 
is close to one. 
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FIG. 4: Electron heating efficiency, a) Schematic representation of ultrafast energy relaxation after photoexcitation, 
through phonon emission, which leads to lattice heating (top), and through carrier-carrier scattering, which leads to carrier 
heating (bottom). In the case of lattice heating, a larger photon energy leads to a larger phase space to scatter to and therefore 
more energy is transferred to the phonon bath, predicting sub-linear scaling of the photovoltage (normalized by absorbed photon 
density) with photon energy. In the case of carrier heating, a larger photon energy leads to a hotter electron distribution (see the 
smeared Fermi-Dirac distributions next to the Dirac cones), predicting linear scaling of photocurrent (normalized by absorbed 
photon density) with photon energy (in the weak heating regime), b) The internal quantum efficiency (IQE), which represents 
the photovoltage normalized by absorbed photon density, as a function of initial electron energy after photoexcitation for 
ambient temperatures To = 40, 100 and 300 K. Error bars are calculated from independent measurement scans and represent 
the 68% confidence interval. The linear scaling through the origin shows that heating dominates the ultrafast energy relaxation 
and therefore that electron heating is efficient. 





















